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Abstract 

We highlight the conceptual issues that arise when one applies 
the quasi-Hermitian framework to analyze scattering from localized 
non-Hermitian potentials, in particular complex square- wells or delta- 
functions. When treated in the framework of conventional quantum 
mechanics, these potentials are generally considered as effective the- 
ories, in which probability is not conserved because of processes that 
have been ignored. However, if they are treated as fundamental theo- 
ries, the Hilbert-space metric must be changed. In order for the newly- 
defined probability to be conserved, it must differ from the standard 
one, even at asymptotically large distances from the scattering centre, 
and the mechanism for this is the non-locality of the new metric, as 
we show in detail in the model of a single complex delta function. 
However, properties of distant bound-state systems, which do not in- 
teract physically with the non-Hermitian scattering potential, should 
not be affected. We analyze a model Hamiltonian that supports this 
contention. 



1 Introduction 

Since the seminal paper by Bender and Boettcherjilj on PT-symmetric but 
non-Hermitian Hamiltonians, subsequent research has gone through various 
stages. First came an exploration of the reality or otherwise of the spectrum 
of non-Hermitian generalizations of well-known soluble models (see Ref. [2] 
for a systematic approach). However, it soon became apparent that some- 
thing additional to the reality of the spectrum was needed if such models 
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were to be viable as quantum theories with a proper probabilistic interpre- 
tation. This is because the natural metric for PT-symmetric models gives 
the overlap of two wave-functions i/j(x) and ip(x) as J dxip*(—x)ip(x), rather 
than the usual J dx(p*(x)ip(x). Since the corresponding norm is not positive 
definite, the theory endowed with this metric does not represent a physical 
framework for quantum mechanics. Instead it turns out to be possible to 
construct [HI IH E] an alternative, positive-definite metric r], which is dynam- 
ically determined by the particular Hamiltonian in question (see Eq. (I12j) 
below). It was also shown [6] that r) provided a similarity transformation 
from the non-Hermitian H to an equivalent Hermitian h, and this equivalent 
Hermitian Hamiltonian was subsequently constructed, frequently in pertur- 
bation theory only, in a variety of models[71 El E]- 

These investigations were all concerned with isolated non-Hermitian sys- 
tems, but more recently some attention has been given to situations where 
a non-Hermitian system interacts with the world of Hermitian quantum 
mechanics. For example, Ref. [10] examined a non-Hermitian analogue of 
the Stern- Gerlach experiment in which the role of the intermediate inho- 
mogeneous magnetic field flipping the spin is taken over by an apparatus 
described by a non-Hermitian Hamiltonian. This type of set-up has been 
further elaborated by Assis and FringfTT] and Guenther al.|12j. although 
the analysis has been challenged by Martin [13] and Mostafazadehp2EJ. In 
another investigation |15] we show in a number of examples that when a 
non-Hermitian Hamiltonian with a real spectrum is coupled to a Hermitian 
system the resulting combined system also has a real spectrum provided that 
the coupling constant does not exceed a critical value. 

It is in this spirit that we reconsider here the conceptual issues involved 
in scattering from a local non-Hermitian potential, first addressed in a some- 
what different context by Znojil[16j, taking as specific examples a pair of 
PT-symmetric complex 5 functions, a finite PT-symmetric complex square 
well[T7l [IH [19], and a single complex 6 function [20]. 

In the next section we review the elementary calculations that show that 
for such potentials the probability, as usually defined, is not conserved. This 
was not seen as a problem by the authors of Refs. [T71 [TU], who presumably 
took the point of view that these potentials were to be treated as effective 
potentials, as is frequently done in nuclear physics, the complex nature of 
the potential representing processes not taken into account, where typically 
the particle can be absorbed rather than merely scattered. 

However, if one takes these potentials seriously as fundamental interac- 
tions, then the whole tenor of the work summarized above is that one should 
change the Hilbert space by adopting the appropriate metric. Of necessity, 
this metric must differ from the standard one not only in the vicinity of the 
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non-Hermitian potential, but also at distances remote from it, in order that 
the newly-defined probability be conserved. 

In Section 3 we reanalyze the case of the single complex delta function in 
the quasi-Hermitian framework. We will see how this gives a mathematical 
resolution of the problem, and it is of some interest to see how it works in 
detail. The essential mechanism is the non-locality of the metric. 

In the final section we address the conceptual issues arising from these cal- 
culations, and in particular emphasize the special status of the plane waves, 
which are substantially affected by the change in the metric, as opposed 
to bound-state wave functions, where a model calculation suggests that the 
changes induced by the introduction of a distant non-Hermitian potential are 
exponentially suppressed. 

2 Some Simple Complex Scattering Poten- 
tials 

2.1 Two delta functions 

The simplest such potential to consider, designed to be PT-symmetric, is a 
pair of imaginary 5 functions: 



V(x) = i\(5(x - a) - S(x + a)) 



(1) 



with a wave-function of the form 




x < —a 



x < —a 



(2) 



Applying the continuity conditions 



(3) 



at each of the interfaces, we readily find that 



1 



D = 



1 + 2ia 2 e 2tka sin 2ka 



(4) 



C = 2iDa(l — a) sin2/ca, 
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where a = A/ {2k). 

Calculating \D\ 2 and \C\ 2 , which in conventional quantum mechanics repre- 
sent the respective probabilities for transmission and reflection we find that 

|n|2 = \ (k) 

1 1 l-Aa 2 (l-a 2 )sin 2 2ka K ' 

which is always greater than 1 for a < 1, as of course is the total probability 

\ r \2 | in|2_ l+4a 2 (l-a) 2 sin 2 2fca 

11 11 l-4a 2 (l-« 2 )sin 2 2fca' W 



2.2 Imaginary Square Well 

A somewhat more complicated potential, still PT-symmetric, is the imagi- 
nary square-well potential [T71 [TBI [19] 

\x\ > a 

-a < x < , (7) 
< x < a 

for which a physical realization, as an effective potential, was given by 
Ruschhaupt et al. [21J. The expressions for the standard transmission and 
reflection coefficients are now considerably more involved |17l [TP"] but Fig. 1, 
which is a graph of the transmission coefficient T for a = A = 1, shows that 
it exceeds unity for a wide range of values of k. 



V{x) 





Figure 1: Standard transmission coefficient of Eq. ([7]) for a = A = 1 
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2.3 Single Complex Delta Function 

Here we consider a single complex delta-function potential of the form 

V = z5(x) } (8) 
where z is complex. The wave-function now takes the simple form 

f e ikx + Ce~ ikx x < , s 

with 

Writing z = 2A(1 + ie) and A; = Xq, it is easy to show that 



2eq 



\C\ 2 +\D\*=[1- , ; ^ ) , (11) 



-1 



which is greater than or less than 1 depending on the sign of e. 

3 Quasi-Hermitian Analysis of the Single Delta 
Function 

The positive-definite Hermitian metric operator r\ = is defined by 

H ] = rjHrj- 1 . (12) 

We are then instructed to calculate matrix elements by including 7]. Thus, for 
example, the expectation value of the quantity represented by the operator 
A in the state \ip) is 

(A)^ = (VMM (13) 

assuming that ip is normalized: (^l^l^) = 1. An operator A is an observable, 
with real eigenvalues, if it, like the Hamiltonian, is quasi- Hermitian, satisfying 
A 1 " = r}Arf l . 

Although the metric operator has been worked out for the complex square- 
well potential by Mostafazadeh[18j the calculation of the corrected transmis- 
sion and reflection coefficients is not particularly transparent, because of 
the complicated nature of the wave function. Instead we will analyze the 
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single complex S function, which, although not PT-symmetric, still has a 
completely real spectrum and therefore falls within the framework of quasi- 
Hermiticity. In this case 77 has again been calculated by Mostafazadeh[20j up 
to third order in the imaginary part of the coefficient of the 5 function. 
In that reference 77 is given as a perturbation series in e, 



V 



X>y r \ ( 14 ) 



r=0 



up to order e 3 . The matrix elements of 7/ ) and rf 1 ' are 

rf°\x,y) = 8{x-y) (15) 

r] [1 \x,y) = ^i\[6(xy)e~ xlx - yl +6(-xy)e- xlx+yl }sgn(y 2 -x 2 ) (16) 

If one is to work with this metric, x is no longer an observable. Instead 
the observable for position is X, defined by 

X = p-'xp, (17) 

where p = Jr\ = e~®l 2 . Thus, indeed, X^ = pxp' 1 = p(pXp~ 1 )p = r/Xr/^ 1 , 
as required. 

The most important property of p is that it is not local (note in particular 
the term in exp[— X\x + y\}). Its effects are felt not only in the vicinity of 
the delta function, but at large distances. Indeed this must be the case if 
unitarity is to be restored within the new framework. 

Of course it has been known for some time that in many cases the metric 
operator is non-local, such as for the potential ix 3 [HE], an d the infinite PT- 
symmetric square-well potential [221 [23]. However, in such cases the particle 
is confined, with an exponentially decreasing wave function so that the range 
of the non-locality is limited. Scattering potentials highlight this feature to 
its full extent because the wave functions of the energy eigenstates do not 
have compact support. 

Let us analyze the scattering problem for the single complex delta func- 
tion in the quasi-Hermitian framework. According to Most af azadeh [2U| |2"2"] 
the relevant wave-function is not ip(x) = (x\if)), but 

= = (x|/#> (18) 

According to this definition, the quasi-Hermitian expectation value of the 
observable X is 

= (^l/T 1 ^ (p- 1 xp)p- 1 \^) (19) 
= = / x\^(x)\ 2 dx 
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The probability density is then taken as g = |\&(x)| 2 , and the total proba- 
bility is conserved[22] in time: d(f g dx)/dt = 0. However, because of the 
non-locality, there does not seem to be a corresponding local conservation 
equation of the form dg/dt + dj/dx = 0. Instead one can show that 

= J dy [**{x)h(x, yMy) - **(y)h(y, x)*(x)] , (20) 

where h is the equivalent Hermitian Hamiltonian h = pHp^ 1 . 
According to Eq. ( fl8|) . we have to calculate 

*(ar) = J dyp{x,y)ij;(y), (21) 

where, to first order, p = r/^ + ^er]^\ and the original wave-function ip(x) 
is given in Eqs. © and (TTUT) . 

The result of this calculation is that for x > 0, and to 0(e), 

e\k 

#>(ar) = De lkx + ^ - — (e~ lkx - (C + D)e ikx ) + 0(e~ Xx ) (22) 

Since we are concerned with asymptotic states, we will neglect the last term. 
One extremely important feature of this equation is that *&>(x) no longer 
represents a pure outgoing wave oc e lkx , but now contains an 0(e) component 
of an incoming wave as well. Thus the physical picture of the scattering is 
completely changed. 

However, in the calculation of probabilities to 0(e) we can neglect this 
term, since it gives rise to a contribution of 0(e 2 ). Doing so, and using the 
simplified notation k = qX introduced earlier, Eq. ( |22|) reduces to 



( x ) = e ikx ^— (1+ , £ . | , (23) 
>y ! q + i\ 2{q + i)' 



so that 



The corresponding calculation for x < 0, again ignoring exponentially 
damped terms, this time oc e Xx , gives 



2(A 2 + k 2 ) J \ 2(A 2 + k 2 ) 

(25) 



+ 1)7 g + z V 2(g + i 
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so that 

\^<(x)\ 2 = { 1 + -t^t) + ~ ( 1+ —5—— —) + interference term, (26) 

V <r + iy g^ + i\ <r + 1 / 

in which the first term represents the incoming flux and the second term 
the reflected flux. Thus we see that each Hermitian flux, for the real delta 
function 2\8(x), is multiplied by the common factor 1 + e/(q 2 + 1), so that 
the newly-defined probability is conserved to this order. 

4 Discussion 

Eqs. (!22l and (1251) provide the resolution of the lack of conservation of prob- 
ability in the normal quantum mechanical framework, which is achieved by 
recognizing that the observable X differs from the usual position variable x, 
even at infinity. Recall that in the calculation to 0(e) we were able to ignore 
the term oc e~ lkx in However, it would have to be taken into account 

in higher orders, and represents a completely different scattering situation 
from the original one, where the wave-function ip(x) represented a wave in- 
coming from the left and being either reflected or transmitted. Instead, the 
corresponding represents a wave with a component also incoming from 
the right. 

However, the scattering problem should be posed in terms of instead 
of ijj(x). Unfortunately this can only be done in a perturbative fashion. That 
is, we have to start with x because that is how the potential is defined. 

The issues we have raised are generic, and not due to the fact that the 
specific potentials we have considered involve 5 functions or sharp discon- 
tinuities in their derivatives. Thus if we consider a general PT-symmetric 
potential V(x) = V*(— x), the integrated continuity equation that results is 

= -2 / lmV(x)\iP(x)\ 2 

—00 J -co 

= / ImV(x)Mx)\ 2 -m-x)\ 2 ), (27) 
Jo 

where j = dip / 'dx — ipdip* /dx) is the standard probability current. 

Barring exceptional cancellations, the right-hand side of this equation van- 
ishes only when ip*(—x) = ip(x), i.e. when the wave-function is itself PT- 
symmetric. However, in the standard scattering set-up whereby a wave enters 
from the left and is either reflected or transmitted ip is definitely not PT- 
symmetric, and j is not conserved. Thus any quasi- Hermitian description 
must necessarily involve a redefinition of x at infinity. 
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The reason the effect of the non-Hermitian potential is felt at large dis- 
tances is that the plane wave functions are completely spread out in position, 
and indeed physically interact with the potential. However, there would be 
a problem if the change in the metric caused by the introduction of the po- 
tential were to appreciably change the properties of particles bound by other 
localized potentials with support far away from it. 

As an example, consider the PT-symmetric potential 

V(x) = -2a5(x) + i\(5(x - L) - 8{x + L)), (28) 

made up of a Hermitian binding potential and two non-Hermitian scattering 
potentials at a distance L, which we will ultimately take to be large. For 
A = the system has a bound state with wave- function ip( x ) ^ e~ a ^ and 
energy a 2 . While the scattering properties of this potential will undoubt- 
edly be affected by the presence of the distant non-Hermitian potentials, we 
would not expect the bound-state energy or wave-function to be appreciably 
changed. 

Writing the wave function as 

Ae KX x < -L 

, Be KX + Ce- KX -L<x<0 

1 = * De~ KX + Ee KX < x < L ' ^ ' 

Fe' KX L < x 



we find 

B 



A 




A 


iX 






2k 6 



2k / 

C = A^- e- 2KL (30) 



D = F[l + lX 



and similarly 



2k 

E = -F — e~ 2KL (31) 
2k 

Finally, applying the continuity conditions at x = we obtain the eigenvalue 
equation for k, which is most conveniently written as 

6 -2kL 

a = K + ? \2a + (K + a) e~ 2KL ] . (32) 

1 + (2/t/A) 2 
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For large L this reduces to 

K o p -2aL 

- ~ 1 - — ff , (33) 

a l + (2a/A) 2 V ; 

showing that indeed the change in the energy is exponentially small. 

From Eqs. (1291) and ( 13TT) we see that although the wave function contains 
an exponentially growing term in < x < L, its coefficient E is suppressed 
relative to D by a factor e~ 2/tL . Thus it only becomes comparable with 
the leading term De~ KX at distances of the order of L, when both terms 
are exponentially small. The wave-function is in fact PT-symmetric, with 
ip(—x) = e t{p il}*(x), where tp is an arbitrary phase. In particular \A\ = \F\, a 
condition equivalent to the eigenvalue equation. 

We have not attempted to calculate rj for this example. It is, of course, 
irrelevant to the calculation of the energy, which is invariant under the sim- 
ilarity transformation induced by rj. It would undoubtedly have the same 
feature of non-locality as that of the previous section, and would necessarily 
induce a large change between ip(x) and ty(x) for the scattering problem, 
where the standard probability current is not conserved. However, as antic- 
ipated, the corresponding changes in the bound-state wave-function would 
be exponentially small for large L because they would involve values of ip at 
large distances from the origin. 

A plausible generalization from this example is that the treatment of 
bound-state problems is essentially unaffected by the introduction of distant 
non-Hermitian scattering potentials, and that it is only when the system 
physically interacts with those potentials that a drastic change in the for- 
malism is required, if indeed the non-Hermitian potentials are regarded as 
fundamental rather than effective. In principle this change provides a consis- 
tent picture of the scattering, although in practice it seems extremely difficult 
to carry out, since r\ can typically be calculated only up to some low order 
in perturbation theory. 
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